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Abstract 

We study transition form factors for radiative and rare semi-leptonic 5-meson decays into 
light pseudoscalar or vector mesons, combining theoretical constraints and phenomeno- 
logical information from Lattice QCD, light-cone sum rules, and dispersive bounds. We 
pay particular attention to form factor parameterisations which are based on the so-called 
series expansion, and study the related systematic uncertainties on a quantitative level. In 
this context, we also provide the NLO corrections to the correlation function between two 
flavour-changing tensor currents, which enters the unitarity constraints for the coefficients 
in the series expansion. 
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1 Introduction 



At the upcoming LHC experiments (notably LHCb), exclusive 5-meson decays will play 
one of the major roles for precision tests of the flavour sector in the Standard Model (SM) 
and its possible New Physics (NP) extensions [TJ[2] . In order to extract information about 
the underlying short-distance flavour transitions, hadronic matrix elements will be required 
as theoretical input, and the precision to which they can be predicted will be essential for 
the success of the flavour program at LHCb. The most prominent example are transition 
form factors (FFs) for 5-meson decays into light mesons, which will be the subject of 
this work. Being non-perturbative hadronic quantities, the theoretical calculation of FFs 
requires techniques such as Lattice QCD (see e.g. |3]-[7]) or QCD sum rules on the light cone 
(LCSR, see e.g. |8]|9], and 10 and references therein) . The two methods are complementary 
with respect to the momentum transfer q 2 between the initial and final-state mesons: In 
Lattice QCD, results are more easily obtainable at high values of q 2 , as discretisation effects 
can only be controlled for small momenta of the final state in units of the Lattice spacing. 
This is in contrast to the LCSR method, which involves an expansion in inverse powers of 
the energy of the light daughter meson that is valid for low values of q 2 . 

Different FF parameterisations, which can be used to interpolate between the results for 
small and large momentum transfer, have been suggested in the literature; a good review 
can be found in j6j. These include simple pole-type parameterisations, like the Becirevic- 



Kaidalov (BK) approach 11 , or variants like the Ball-Zwicky (BZ) parametrisation [8J. 
Another representation is found from the Omnes solution to the dispersion relation, see 
the discussion in [l2lfl3]. In this paper, we will make use of the so-called Series Expansion 



(SE), which was advocated in Refs. 14 -18 . Here, one can make use of dispersive bounds to 
obtain additional theoretical constraints on the expansion coefficients. A simplified version 
of the SE (SSE) was recently suggested in [L3^ . The aim of this work is to use the SE/SSE 
to describe the transition FFs on the basis of recent Lattice and LCSR results, including 
a detailed analysis of systematic errors. We will focus in particular on the FFs entering 
B — > Vj, B — > L£ + £~ , B — > Lvv decays, where L = P, V is a light vector or pseudoscalar 
meson. We will give numerical results for B — >■ p, K, K* and B s — > <fi transitions, which 
are particularly interesting with respect to NP studies, see e.g. [H[2| [l9}j23] . (Detailed 
phenomeno logical studies for FFs relevant for the determination of the CKM elements \ V u b\ 
and \Vcb\ from semi-leptonic B decays can be found in the recent literature [3}]5| [l3|[24] .) 
For the discussion of dispersive bounds for tensor FFs, we will include the result of a precise 
calculation of the tensor current two-point correlator at NLO in the QCD coupling constant, 
including the leading non-perturbative corrections from quark and gluon condensates. 

Our paper is organized as follows. In Section [2j we provide convenient definitions for the 
.B-meson FFs and introduce the idea of the SE/SSE. In Section |3j we review the derivation 
of dispersive bounds from current-correlation functions and summarize the results for the 
profile functions obtained from the operator-product expansion. We apply our formalism 
to B K, B -»■ p, B -»■ K* and B s -»■ <\> FFs, by fitting the (truncated) SE/SSE to 



theoretical "data" from Lattice QCD and/or LCSRs in Section 4.3 Our conclusions are 



presented in Section [5j and some technical details are given in the Appendix. 
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2 Form Factors 



In this section, we will first provide the definitions for the various I?-meson FFs in question, 
fix the notation to be used in the subsequent discussion, and introduce the SE/SSE. 

2.1 Definition of Form Factors and Helicity Amplitudes 

The hadronic matrix elements for transitions between a pseudoscalar B meson and a generic 
(light) pseudoscalar meson are usually written in terms of three FFs fo(q 2 ), f+(q 2 ) and 
frig 2 ), which depend on the momentum transfer q 2 = (p — k) 2 , 

(P(k)\q lfl b\B(p)) = (p, + ^ ~ 2 ^ Utf) + % fo(q 2 ) , 

(P(k)\qa, u q%\B(p)} = % - (q 2 (p + - (m| - m%) gj f T (q 2 ) • (1) 

m B + Tflp 

At zero momentum transfer, the additional relation /+(0) = /o(0) holds. 

Similarly, the matrix elements for the transition between a S-meson and a generic vector 
mesoiiQcan be written in terms of FFs V(q 2 ), Ao_3(g 2 ), Ti_3(g 2 ), which are conventionally 
defined as 

(V(k,e)\q^b\B(p)) = ie^ pa e*»(k)p<>k° 21 ( ' /_) 



(V(k,e)\q W b\B(p)) = -s;(k)(m B + m v )A 1 (q 2 ) + (p + k),(e*(k)-q)^f' 



+q^e*(k)-q) 2 ^(A 3 (q 2 )-A (q 2 )) , (2) 

where Ao(0) = A^O). For transitions involving a tensor current, the matrix elements are 
characterised by the tensor FFs, 

(V(k,e)\qa^q v b\B(p)) = i w s* v tfk° 2T x {q 2 ) , 

(V(k,e)\qa^ 5 b\B(p)) = T 2 (q 2 ) (e^k) (m| - m\) - (e*(k) ■ q) (p + 

-,2 



+T 3 (q 2 )(e*(k) ■ q) (q, - _ 2 ~ _ 2 (2p - q), ) , (3) 



m B -my 



where Ti(0) = T 2 (0). The equations of motion for the quarks imply the additional con- 
straint 



2 m B + m v . , ox m B — mv A , 2 \ 
and therefore the B — y V transitions are characterized by seven independent FFs. 



Our phase convention for the vector state differs by a relative factor of i from the convention that is 
used, for instance, in [9j. 
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The above constitute the standard definitions for the FFs widely used in the literature. 
However, for this work, we find it convenient to work with certain linear combinations of 
these, dubbed helicity amplitudes in Ref. 16 . This is primarily to diagonalize the unitarity 
relations, which shall be used to derive the dispersive bounds on certain FF parameteri- 
sations. The helicity amplitudes also have definite spin-parity quantum numbers, which is 
useful when considering the contribution of excited states. In addition, they have simple 
relations to the universal FFs, appearing in the heavy-quark and/or large-energy limit (see 
appendix [Bj , and lead to simple expressions for the observables in B — > L£ + £~ decays in 
the naive factorization approximation. To put the contributions to the various correlation 
functions entering the dispersive bounds on an equal footing, we also choose a particular 
normalization convention and define new B — > P vector FFs via 

AvAl 2 ) = (P(k)\qj»b\B(p)) • (5) 

Here 

A = ((m B - m P f - q 2 ) ((m B + m P ) 2 - q 2 ) = (t_ - q 2 ){t + - q 2 ) (6) 

is a standard kinematic function, and e*J l {q) are transverse (a — ±), longitudinal (a = 0) 
or time-like (a = t) polarization vectors as defined in (60) in the appendix. This implies 

A Vfi (q 2 ) = f + (q 2 ) , AvM 2 ) = "^f^ /o(g 2 ) , (7) 
while the transverse projections vanish. Similarly, for the B — >■ P tensor FF, we define 

A T Aq 2 ) = H)Jle*S(q) (P(k)\qa^ b\B(p)) ■ (8) 
Here, the only non-zero FF ia^] 

A T ,o(Q 2 ) = fr(q 2 ) ■ (9) 

A similar analysis for the B — > V vector and axial-vector FFs yields 

BvAq 2 ) = M (v(M*))|s7m(i ~rW0°)> ( 10 ) 

e(fc) 



2 The newly defined tensor FF AT,o(q 2 ) vanishes as yV, which might look somewhat artificial at first 
glance. However, the tensor current does not contribute to physical processes at q 2 — anyway. 
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with 



2\ 



(m B + m v ) 2 (m% - m v - q 2 ) A x (q 2 ) - A A 2 (q 2 ) 



2m v \f\ (m B + m v ) 
By t {q 2 ) = A (q 2 ) , 

Bv,i(q = ^ — = — , v (i 

m B + m v 



B v ,- 








B v> - 


- + By + 


V2 



Bvfl{<?) = — ~ v ~ v^^^ A^). (11) 

v A 



Finally, the B — > V matrix elements with tensor currents are projected on 

BtM) £ T(ti (V(k 7 e(k))\qa, a q a (l+j 5 )b\B(p)) (12) 



giving rise to the FFs 

n (a 2 \ - Vq 2 ( m % + 3m v - q 2 ) rp (n 2s Vq 2 ^ T , 2] 

B T Aq 2 ) = - Bv '-^ Bv ' + =V2 Tl (q 2 ), 

BtM 2 ) = - ^ ( ^-<W )- d3) 

2.2 Series Expansion 
2.2.1 Resonances 

An important factor in determining the shape of the FF is the presence of low-lying res- 
onances with appropriate quantum numbers and mass in the range £_ < < t + . 
Common to most parameterisations is the inclusion of the low-lying resonance by a sim- 
ple pole. The various descriptions differ in the modelling of the continuous part. In the 
following, we use the abbreviation P(q 2 ) = 1 — q 2 jm\. If multiple resonances are present 
in the given region, then P(q 2 ) should be a product of such poles, and if no resonances 
are present then P(q 2 ) = 1. A summary of the relevant resonance masses is provided in 
Tabled! 



2.2.2 Series Expansion(SE) 



The starting point is to extend the 
^ to analytic functions 



The SE has its origin in dispersive relations 14 16 18 
FFs defined in the physical range (from q 2 
throughout the complex t = q 
production of real BP/BV pairs at q 2 > t + = (m B + mi) 2 - If low-lying resonances are 



to t- = \m B — mLj 
2 plane, except for along the branch cut at the threshold for 
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Table 1: Summary of the masses of low-lying B^ and B s resonances, using PDG values 
25 and/or theoretical estimates from heavy-quark/chiral symmetry [26). Notice that 
the mass values for (0 + , 1 + ) predicted in 26 have not been confirmed experimentally, yet. 



Instead the PDG quotes "effective" resonances Bj (5698) and -B*j(5853) with undetermined 
spin/parity. 



Transition 


J p 


Mass (GeV) 


J p 


Mass (GeV) 


Ref. 


b -> d 


0" 


5.28 


\- 


5.33 






0+ 


5.63 


1+ 


5.68 






1+ 


5.72 


2+ 


5.75 


25 


b -> s 


0" 


5.37 


\- 


5.42 






0+ 


5.72 


1+ 


5.77 






1+ 


5.83 


2+ 


5.84 





present below t + , they are accounted for by the so called Blaschke factor B(t), see below. 
Complex analysis can then be used to map the cut t-plane onto the unit disc in terms 
of the coordinate z(t). The variable z(t) is found to be an excellent expansion parameter 
for the FFs. Furthermore, with an appropriately chosen normalization function 0/(t), one 
obtains simple dispersive bounds on the coefficients of the SE, see below. We will discuss 
the calculation of the functions </>/(£) as well as the derivation of the dispersive bounds in 
greater detail in Section [3] The Series Expansion (SE) then corresponds to the following 
FF parametrisation, 

'M = 5(i)b> !>**<«>■ (14) 



with 

z(t) = z(t, to) = ^f=f^f=f • («> 

Here < to < t- is a free parameter which can be optimised to reduce the maximum value 
of \z(t)\ in the physical FF range, 



*oL =U[l-Jl-^). (16) 



opt. t ^ y t 

We will later see that with the optimised value for to, the FFs can be well described by a 
SE which is truncated after the second term proportional to z(t). Other values of to (e.g. 
to = 0) are still allowed but sometimes require to go to higher order in the SE. 

As a crucial property, the function z(t) satisfies \z(t)\ = 1 in the pair-production region, 
t > t + . The Blaschke factor is thus chosen as B(t) = z(t, m^). As for P(q 2 ) defined above, 
if multiple resonances are present, then B(t) is a product of the corresponding Blaschke 



factors. Further discussion about the physical basis for the SE is found in Ref. 16 
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2.2.3 Simplified Series Expansion (SSE) 



Another form of the Series Expansion method can also be considered. Instead of the 
Blaschke factor B(t), one can use a simple pole P(q 2 ) to account for low-lying resonances. 



This idea was proposed in Ref. 13 , yielding 



It was found that the dispersive bounds can still be imposed on the coefficients at of the 
SSE. We will discuss this and other issues concerning the validity of the simplifications in 
the following section. 



3 Dispersive Bounds 

The FFs describe the process B — > L with L = P, V in the decay region < q 2 < 
£_ = (mg- m^) 2 . Using crossing symmetry, they can also describe the process in the pair- 
production region q 2 > (mB+m L ) 2 . This can be exploited to obtain a bound on parameters 



describing the FFs. A detailed derivation of this bound can be found in Refs. 15,17 . Here 
we provide a brief outline of the argument, in order to introduce our notation and to extend 
the method to tensor FFs. 

The crucial observation of the idea of dispersive bounds (as it is for QCD sum rules) is 
the possibility to evaluate the correlator of two flavour-changing currents, 

K(l 2 ) =i J d'xe^(0\T^(x)jl x (0)\0) , (18) 

either by an operator product expansion (OPE) or by unitarity considerations. Here the 
relevant currents are defined a^] 



% = rtj>, j'r A = ^(i-7 5 )& 



3% = > lT AT = + l 5 )b . (19) 

Furthermore, we introduce longitudinal and transverse helicity projectors, 

W) = > = ^zi(^f- r) , (20) 

which allow us to rewrite the correlation functions in terms of Lorentz scalars, 

Tlf(q 2 )=Pr(q 2 )K(q 2 ), (I = L,T). (21) 



3 In phenomenological applications, we are only interested in the currents j'J +At • The connection to 
correlators with genuine tensor currents — qa^ u q is given in Appendix |T)| 
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As Uf(q 2 ) is an analytic function, it satisfies the subtracted dispersion relation, 

(22) 



n! (io 2?1 



(t-g 



2 \n+l 



<? 2 =0 





A'/ 



where the number of subtractions n is chosen to render the resulting function Xi( n ) finite. 
3.1 Hadronic representation of the Correlator 

Unitarity allows us to express ImlTf (g 2 ) as the positive definite sum over all hadronic 
states T with allowed quantum numbers: 

ImlW) = \ W dp T (2tt) 4 5\q - p r ) If (0| j* \T) (T\ £* |0> . (23) 
r ^ 

where pr is the total momentum of the final state, and dpr contains the appropriate phase- 
space weighting. For a particular choice of intermediate state, T = BL, we define 

Imn^(? 2 ) = rj J dp BL If <0|# \BL) (BL\jI x |0> , (24) 

where rj is an isospin-degeneracy factor for a given channel, and we relegate the contribution 
from phase space to the function 

1 f d 3 p B d 3 p L 4 



dpBL = ^ J 2E- B 2E- L nq ~ PB - pL) - (25) 
Clearly, this results in the inequality 

lmU x BL (t) < Imnf (*) . (26) 

Now, by extending the FFs to analytic functions throughout the t-plane, except for along 
the branch cut at the threshold for production of real BL pairs, one can use crossing sym- 
metry to relate the matrix elements (0| j x \BL) to (B\ j x \L). The latter can be rewritten 
in terms of FFs, as defined in Section [2] As stated earlier, we use helicity-based linear 
combinations of the traditional FFs, such that all production amplitudes 'diagonalize': 



(P\rt\B)(B\jt v \P) = ^lAol 2 

A 



(P\fi\B)(B\ji v \P) = - 2 \Av, t \\ 
Pf (P\fi\B)(B\jt T \P) = ^\A Tfi \\ (27) 
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for B decays into pseudoscalars, and 



Pf(V\jV- A \B)(B\j}; v - A \V) 



pr(v\j:- A \B)(B\jtv- A \v) 



Pf (v\£ +AT b\B)(B\ 3 t> T+A n\v) 



3q~ 2 



El*. 



2 

1 1 j 



i=0 



A \n i 2 
g 2 



E i e « 

i=Q 



We can now express Im Hf BL in compact form 
lmUf BL = r] J dp BL j t \A* 



X I 4x12^ g A 3 / 2 ■ x , 2 
" 48vr t 2 1 7 1 



where the \Af\ can be read off (27 28) 

VI 2 i a \2 



A 



\>Av,o\ 



A 



V|2 



3U 



v,t| 



for decays into pseudoscalars, and 



El"-. 



2 

1 1 ) 



K- A r = 3|Bv;*l 



i=0 



= q 2 \A T ,o\ 



\A T T +A -\ 2 = q *Y,\®T. 

i=0 



i2 
i I > 



for decays into vector mesons. 



(2* 



(29) 



(30) 



(31) 



3.2 OPE for the Correlator 



Alternatively, we can examine the correlator (18), using an OPE for the T-ordered product 
of currents in the limit q 2 = <C t + . The standard expansion takes the form 27 29 



/CO 
dxe^PrT{jX(x)ji x (0)} = J2 C? ik (q)<D, 
k=i 



(32) 



where Cf n (q) are Wilson coefficients for a given current X and projector /, and O n are 
local gauge-invariant operators, consisting of quark and gluon fields. Here, the operators 
are ordered by increasing dimension k. We can use the above, to express the correlator, 



(33) 



k=l 



Besides the identity operator, whose Wilson coefficient contains the purely perturbative 
contribution to the correlator, we will specifically consider the first few operators related to 
the non-perturbative contribution from the quark condensate (m q qq) , the gluon condensate 



s G 2 ), and the mixed condensate (g s q(a • G) q). We will elaborate on our calculation of 

the Wilson coefficients, Cf k (q 2 ), later. Specifically, we must calculate the Wilson coeffi- 
cients entering the functions xf{ n ) i n (22). 



7T 
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3.3 Bounds on coefficients in the SE 



Using (26), we find 



1 

7T 



dt 



imn^ L (i) 



(t - g 2 ) 



2 \n+l 



1 

7T 



77 A 3 / 2 (t) 
48vr t n + 3 



^?(t)r<xf(n), 



(34) 



*+ 



where xf = X/ope * s calculated from (33). Mapping the pair-production region t > t H 
onto the unit circle |;z(t)| = 1, this inequality could be written in the form 



1 

2vri 



dz 

z 



^Af\ 2 [ 



z) < 1 



1 

7T 



* — *n 



*-t. 



°'-Mf| 2 (0<l, (35) 



where the function (t)| 2 can be obtained by comparing (35) and (34), and using \{t) 
(t+-t)(t_-t), 



I0?(*)P 



(t - *+) 2 (t-t-) 3/2 t-t 



487rxf (n) (t+ - to) 1/2 t n+2 



(36) 



The isospin-degeneracy factor 77 takes the values 3/2, 2 and 1 for B — > p, B — > and 
B s — > <f) respectively. We may now generically write the helicity-based FFs Af(t) as 



Af(t) 



B(t)<Pf(t) 



k=0 



(37) 



with real coefficients c^, and a Blaschke factor B{t) = Yl i z(t, rn 2 R ), representing poles due 
to sub-threshold resonances of masses m^, and satisfying |-B(t)| = 1 in the pair-production 
region. The additional factors (y/—z(t, 0)) m and z(t, t_)) 1 have been added to take into 
account the unconventional normalisation of our FF functions through factors of ^fq 2 and 
a/A (e.g. m — 1 for At,o, and I — — 1 for ^4.y it , cf. above) Q The function <f)f(t) has to be 
constructed in such a way that its absolute value satisfies Eq. (36), while (|37|) retains the 



analytical properties of the FF. This can easily be achieved by replacing potential poles 
and cuts in a/10;^ (t)| 2 , by making replacements of the form 



— > 



-z(t,X) 



(3* 



t-X t-X ' 

which is allowed as \z(t, X)\ = 1 in the pair-production region. This results in (see also (30 



0fW 



77 



(t-t+) fz{t,0) 



487T X f(n) (^-to) 1 / 4 



-t 



(3+n)/2 



z(Mo) 
t -t 



-1/2 



t- ~ t 



-3/4 



(39) 



4 These factors could also be considered as part of the Blaschke factor. Note that under a change of 
normalisation convention for the FFs, both, the so-constructed Blaschke factor as well as the function 0(£) 
have to be modified, while the coefficients of the SE remain the same. 
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Inserting the parametrisation (37) into (35), and using \z(t, to) I = \ z (tj m %) \ — l 2 ^, 0)| = 1, 
the integration dz/z = dip along the unit circle is trivial, yielding the desired bound on the 
coefficients 



XX <i. 



(40) 



fc=0 



For decays into vector mesons, using an analogous parametrisation as (37) for each 
individual FF contribution in (31), one obtains a bound on the sum of the corresponding 
coefficients. As an example, let us consider A^~ A (t) : where we parameterise 



B Vfi (t) 
Bv,i{t) 
BvAt) 



B{t)^J7)<p v T - A {t) 



K-1 
fc=0 



(v,o) z k 



B{t)4> v T - A {t) k M 



K-1 



B(t)v^nj$- A (t) ,„ 

resulting in the dispersive bound 



(V,2) z k 



(41) 



K-1 



E ((^ y,0) ) 2 + i^f + (/^ 2) ) 2 ) < i • 



(42) 



k=0 



3.4 The coefficients xf( n ) 

In Table |2] we summarize the numerical result of our calculation of the various coefficients 
xf{ n )i which enter the functions <f)f{t) in the SE. We quote individual numbers for the 
perturbative LO and NLO results, as well as from the condensate contributions, for two 
different values of light-quark masses, m q = ma and m q = m s . Also the number of sub- 
tractions is indicated. Details of the calculation as well as analytical formulas can be found 
in Appendix [Cj As can be observed from Table |2j the NLO perturbative corrections are 
essential for a reliable estimate for the coefficients xf{ n )i while the quark condensate gives 
only small contributions, and the gluon condensate and the mixed quark-gluon condensate 
are negligible. 



4 Form Factor Fits to Theoretical Data 
4.1 Theory Input from Lattice and LCSR 

There are two main methods to obtain theoretical predictions on 5-meson decay FFs: 
Light-cone sum rules (LCSRs) and Lattice QCD. As mentioned in the introduction, these 
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Table 2: Summary of OPE results for the coefficients xf 
values have been used (27 31 33]: [i 



a s = 0.2185, (dd) 
(1.4 GeV) 2 (qq). 



m b 



(278 MeV) 3 , (ss) 



n). The following parameter 
= 4.2 GeV, m d = 4.8 MeV, m s = 104 MeV, 
0.8 (dd), (^G 2 ) = 0.038 GeV 4 , (qGq) = 



q 


Correlator 


Subtractions 


LO 


NLO 


(??) 




(qGq) 






100 x mlx s 


2 


1.265 


0.589 


0.029 


0.001 


-0.003 


1.88 




100 x m 2 h x P 


2 


1.268 


0.590 


0.029 


0.001 


-0.003 


1.88 




100 x xl 


1 


1.262 


0.211 


0.029 


0.001 


-0.003 


1.50 


d 


100 x xi 


1 


1.271 


0.205 


0.029 


0.001 


-0.003 


1.50 




100 x m\x v T 


2 


0.951 


0.236 


-0.029 


-0.001 


0.007 


1.16 




100 x mlxr 


2 


0.948 


0.237 


-0.029 


-0.001 


0.007 


1.16 




100 x mlxr 


3 


2.539 


0.579 


-0.029 


-0.000 


0.008 


3.10 




100 x itiIxt T 


3 


2.527 


0.586 


-0.029 


-0.001 


0.008 


3.09 




100 x mftX S 


2 


1.233 


0.571 


0.024 


0.001 


-0.003 


1.83 




100 x m\x P 


2 


1.296 


0.608 


0.022 


0.001 


-0.003 


1.93 




100 x x ^ 


1 


1.172 


0.229 


0.023 


0.000 


-0.003 


1.42 


s 


ioo x xi 


1 


1.361 


0.187 


0.023 


0.002 


-0.003 


1.57 




100 x m\x v T 


2 


0.980 


0.237 


-0.022 


0.000 


0.005 


1.20 




100 x m^XT 


2 


0.916 


0.238 


-0.024 


-0.002 


0.006 


1.13 




100 x trIxt 


3 


2.652 


0.569 


-0.023 


0.001 


0.006 


3.21 




100 x mlxT T 


3 


2.404 


0.603 


-0.024 


-0.002 


0.007 


2.99 
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techniques are largely complementary, as they perform best in different regimes of mo- 
mentum transfer q 2 . It is worth mentioning that certain decays, e.g. decays to unstable 
hadrons, are more challenging in Lattice QCD, and in some cases only quenched results 
exist for a subset of the FFs. On the other hand, LCSRs provide results for all decay 
channels considered in this work, including the complete set of seven FFs for B — > K* and 
B s — > (f>, which so far have not been fully addressed by Lattice QCD. However, as LCSR 
results are only valid in the low-g 2 regime, in these cases further theoretical information is 
needed to extrapolate to large values of q 2 , as will be discussed in the following subsections. 

In our analysis, we will use the LCSR predictions from Refs. |8] and j9], taking 3(4) 
points at low values of q 2 as input, see Table [3} The quoted errors are extrapolated from 
the value quoted for q 2 = in the references specified in the Table. Lattice data is available 
for B — >■ p and B — )■ K decays, and is as shown in Table [4|p1 For those data points which 
have an asymmetric statistical or systematic error, in order to perform the fit, we take the 
FF to be the central value in this statistical or systematic range, and take half the range 



to be the statistical or systematic error 34 



For B — )• p and B — )■ K decays, we use LCSR and Lattice data to interpolate between 
the low and high-g 2 region. The result can be compared to the case where we extrapolate 
to the high-g 2 region only on the basis of LCSR predictions. This procedure will gives us 
an idea about the confidence in the extrapolations for those cases where Lattice data is 
lacking. 

4.2 Parameterisation of FFs as Series Expansion 

For those channels where where Lattice data is not available, it is essential to employ a 
FF parameterisation that takes into account the characteristic features of the FF shape as 
determined from the analyticity and unitarity consideration above. For every considered 
FF, we will therefore define a parameterisation based on the SE, 



K-l 



(V,0) k 

a k z 



k=0 



K-l 



AvAt) Bit) yfifcij) m) ~ 



(V,t) k 



k=0 



K-l 



A ^-Bm¥)u ' ( } 

5 We are very grateful to Sara Collins of the QCDSF collaboration for providing us with specific values 
for B -> K. 
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Table 3: Overview of LCSR points used, transformed to the helicity amplitude basis. 



FF 



LCSR/g 2 (GeV 2 ) 



Ref. 



6 



9 



12 



Table 3, 



Av,o 



0.40 ± 0.05 
0.40 ± 0.05 
0.13 ±0.01 



0.48 ± 0.06 
0.51 ±0.06 
0.22 ±0.02 



0.59 ±0.07 
0.65 ±0.08 
0.34 ±0.03 



6 



9 



12 



Table 8, 



B v ,o 
B v ,i 

B V ,2 

B v ,t 
Bt,o 
B t ,i 

Bt,2 



0.37 ±0.12 
0.16 ±0.01 
0.16 ±0.02 
0.37 ±0.04 
0.17 ±0.35 
0.45 ± 0.04 
0.46 ± 0.04 



0.46 ±0.13 
0.27 ±0.02 
0.29 ±0.03 
0.46 ± 0.04 
0.3 ±0.26 
0.55 ±0.05 
0.58 ±0.05 



0.60 ±0.14 
0.41 ±0.04 
0.46 ±0.04 
0.58 ±0.06 
0.47 ±0.23 
0.69 ±0.06 
0.76 ±0.07 



0.71 ±0.22 
0.9 ±0.08 
1.0 ±0.1 



9 



12 



Table 8, 9 



B v ,o 
B v ,i 

B V ,2 

B v ,t 
Bt,o 
Bt,i 

Bt,2 



0.45 ±0.13 
0.19 ±0.02 
0.20 ±0.02 
0.44 ± 0.04 
0.23 ±0.36 
0.59 ±0.06 
0.61 ±0.06 



0.56 ±0.13 
0.32 ±0.03 
0.35 ±0.03 
0.54 ±0.05 
0.39 ±0.27 
0.72 ±0.07 
0.77 ±0.07 



0.73 ±0.15 
0.49 ±0.04 
0.57 ±0.06 
0.67 ±0.06 
0.60 ±0.24 
0.89 ±0.08 
1.0 ±0.1 



0.90 ±0.22 
1.2 ±0.1 
1.4 ±0.1 



6 



9 



12 



Table 8, 



B v ,o 
B v ,i 

B V ,2 

B v ,t 
Bt,o 
B t ,i 

Bt,2 



0.55 ±0.12 
0.2 ±0.02 
0.21 ±0.02 
0.56 ±0.04 
0.26 ±0.39 
0.59 ±0.06 
0.61 ±0.06 



0.68 ±0.13 
0.34 ±0.03 
0.38 ±0.04 
0.68 ±0.05 
0.44 ±0.29 
0.72 ±0.07 
0.77 ±0.07 



0.85 ±0.14 
0.52 ±0.04 
0.62 ±0.06 
0.85 ±0.06 
0.67 ±0.26 
0.89 ±0.08 
1.0 ±0.1 



1.0 ±0.3 
1.2 ±0.1 
1.4 ±0.1 
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Table 4: Overview of Lattice points used, transformed to the helicity amplitude basis. 
Note that specific values for B — » p are as in Table 2 of Ref. 12 . 



q 2 (GeV 2 ) 



FF 



Ref. 



14.5 
15.6 
16.7 
17.9 
19. 
20.1 
21.3 
22.4 



0.94 ±0.19 

1.1 ±0.2 

1.2 ±0.2 
1.4 ±0.3 
1.6 ±0.3 
1.9 ±0.4 

2.3 ±0.4 
2.9 ±0.6 



A v ,t 

1.1 ±0.2 
1.3 ±0.3 
1.5 ±0.3 
1.8 ±0.3 

2.3 ±0.4 
3. ±0.6 

4.4 ±0.8 
8.7 ±1.7 



A 



T,0 



QCDSF 



B 



r.o 



B 



B 



T.2 



UKQCD 



12.7 
13. 
13.5 
14. 
14.5 
15. 
15.5 
16. 
16.5 
17.1 
17.6 
18.2 



0.64 ±0.78 
0.71 ±0.72 
0.8 ±0.66 
0.9 ±0.62 

1.0 ±0.6 

1.1 ±0.6 

1.3 ±0.7 

1.4 ±0.8 
1.6 ±0.9 
1.8 ± 1.2 
2.1 ± 1.5 

2.5 ±2. 



0.34 ±0.27 
0.39 ±0.25 
0.48 ± 0.22 
0.58 ±0.19 
0.68 ±0.16 
0.78 ±0.15 
0.89 ±0.15 

1.0 ±0.2 
1.2 ±0.3 
1.4 ±0.4 
1.7±0.6 

2.1 ±0.9 



0.9 ±0.18 
0.96 ±0.18 

1.1 ±0.2 

1.2 ±0.2 

1.3 ±0.2 

1.4 ±0.2 

1.6 ±0.2 
1.8 ±0.2 
2.1 ±0.2 
2.4 ±0.2 

2.7 ±0.3 
3.3 ±0.3 
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and 



K-l 



(v,t) z k 



K-l 



n (t\ - r (t,o) „ k 



s(t)v^(tX)4 +AT (t) ^ 



fc 2 



^ 2(t) = »m ATTTT^ T^ E ^f > ( 44 ) 



and By_o-2 already given in (41). Here we have used our FF convention defined in 
fl5|8|10fi~2l ) and explicitly quoted the pre-factors, necessary to obtain the correct analyti- 
cal behaviour of our FFsQ In our fits below, we will find that in general the SE can be 
truncated after the first two terms, i.e. the parameter K can be set to 2. We should point 
out, however, that this does not necessarily imply that higher-order terms in the SE are 
negligible: Although \z\ 2 -C 1 in the semi-leptonic region, one may still have [ac^ z 2 \ ~ \a±z\ 
if the coefficients satisfy «i C «2 ^ I- From the theoretical point of view, this reflects 
an irreducible source of uncertainty, which we will discuss in some detail for some spe- 



cific examples in Sec. 4.4 From the practical point of view, we consider the truncated 
SE as a reasonable parameterisation which is easy to implement (and easy to refine) in 
phenomenological studies. 



For simplicity, we will not explicitly implement the theoretical relations ( 64|65 ), that 



some of the FFs fulfil at q 2 = 0, into the fit, because they are automatically satisfied by the 
rather precise input from LCSR at this point. However, the helicity-based FF definition 
further implies a relation between the FFs By a and By 2, and similarly between Bto an d 



Bt,2, see (66) in the appendix, which we will implement as an additional constraints on 
the corresponding coefficients in the SE. From the above parameterisations, the SSE is 
obtained by the replacements 

(j>f(t)^l, B(t)^P{t), y/-z(t,0) V?/^, y/z(t,t-)^VX/m 2 B , (45) 
with new coefficients and (3k- 



6 In Refs. [18] and 13 17 , the predictions from perturbative QCD for the scaling of the FFs at large 
values of q 2 have been used as an additional constraint on the shape of the FFs. We have found that these 
constraints do not influence the FF fits in the decay region significantly. As the asymptotic behaviour of 
exclusive observables in QCD is still a matter of controversy, we therefore find it safer and simpler not to 
include these constraints in our parameterisation. 
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4.2.1 Unitarity constraints 

For the SE parameterisation, the unitarity constraints take the form 

K-l K-l 



( a kf < 1 for A v,o and A T ,o, 3 ^ {a k f < 1 for A Vjt , (46) 



and 



k=0 k=0 



K-l 



fc=0 

£{(/?f 0) ) 2 + (/^P) 2 + (/3P) 2 } < 1 for Bk,o, and S^, 

k=0 

^{(^y + (/3f X) ) 2 + (/?f 2) ) 2 } < 1 for Br,o, Sr,i, and S T , 2 . (47) 

fc=0 

For the SSE parameterisation, imposing the unitarity bound is more complicated, as 
shown in Ref. 13 . We repeat the derivation of this bound in order to define notation used 
later. One first compares the SE and SSE parameterisations: 

K-l K-l 

o k z k =k{z)Y,a k z k (48) 

k=0 k=0 

One can simply obtain A(z) by combining the prefactors from the SE expansion with the 
prefactors from the SSE expansion, and expressing the result as a function of z(t, t Q ). Since 
z is a small parameter, we can expand A(z) in powers of z: 

A(z) = J2(kZ k . (49) 

k 

We can therefore obtain a relation between the coefficients a k and a k , 

min[K—l,i] 

a t = Yl d- k a k) 0<i<K-l, (50) 

fc=0 

which results in bounds of the type 

K-l 

^ C jk atj a k <l, (51) 

j,k=0 

where 

K— 1— max[j',fc] 

(52) 

is a positive definite matrix. 
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4.3 Fitting prescription 

We perform a fit to the LCSR data, as well as, where possible, a combined fit to the LCSR 
and Lattice data, by minimising a x 2 function defined by 

X\6) = (*l - F(U, 0)) [V-% [F 3 - F(t s , )) , (53) 

where 9 contains the parameters of a given FF parameterisation, F t are the FF values from 
LCSR/Lattice at given points % and Vij are elements of the covariance matrix as defined 
below. 

As explained above, we are going to investigate parameterisations based on two variants 
of the SE, where the parameters will be subject to additional constraints derived from 
dispersive bounds on the FFs. 



In the conventional series expansion (SE), we use (14), and truncate the series after 
the first 2 terms, such that 



-» v — -\ ! 

9 = {a , ai} , a i < 1 



The simplified series expansion (SSE) uses (17), with 



9 = {a ,ai}, dj OLjOLj < 1 , 

i,j=0 



where the matrix CL- is defined in (52). 



In constructing the covariance matrix, when we do a combined fit to LCSR and Lattice 
data, we assume the matrix to be block diagonal with independent blocks for Lattice and 
LCSR, equivalent to \ 2 = Xlcsr + xLt» where 

XlcsrW = (F-F(tJ)) [Vj^^Fs-Fte,?)) , (54) 



and 



XLM = (Fi - F(t u 9)) [V^l] . . (Fj - F(t v 9)) , (55) 

We consider the statistical and systematic contributions to the Lattice errors sepa- 
rately. Where results were not available in the literature, we received the breakdown by 
private communication with the authors. In obtaining the covariance matrix, we make the 
following conservative assumptions: 



Statistical errors of Lattice data are 50% correlated [4^7]. 
Systematic errors of Lattice data are 100% correlated [4l[7j 
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Errors of LCSR data are due to parametric as well as to systematic uncertainties from 
different sources. In order to provide a concrete number for the \ 2 value characeter- 
izing the quality of the fit, we have estimated the errors of LCSR data at different 
values U to be 75% correlated^] 



This prescription leads to a covariance matrix V 1 ^ = covft 1 ,^], containing 

and 



1 3 

4 3 4 

Kit = iaVay + ioV + eV 



vy 

LCSR 



(56) 
(57) 



where cxj are the statistical errors, Si are the systematic errors for the Lattice data, and K{ 
are the errors for the LCSR predictions. 

Minimising x 2 ($) then yields the best fit parameters 8*, as well as the covariance matrix 
of the fit, Uij = cov[#j, 9$], 

1 d\\8) 



v 2 dBidBi 



{hi 



from which we calculate the error associated to the fitted FF function: 



AF(t,6*) 



dF(t,9) 
88,. 



U, 



08, 



(59) 



4.4 Results 

Having established the fitting procedure, we consider FFs for the decays B — )■ p, B — > 
K, B —7- K* and B s — > <p. We concentrate on radiative and rare semi-leptonic decays, 
as previously the dispersive bounds had not been calculated for the tensor current, so 
could not be applied to these decays. The phenomenological motivations for studying the 
chosen decays are as follows. First, they involve flavour changing neutral currents via 
e.g. electroweak penguins, so they are particularly sensitive to new physics. Secondly, the 
di-lepton signature can easily be detected at the LHC, and the three-body or four-body 
final state (for subsequent decays K* — > Kir and — > KK) involves many promising 



observables related to various angular distributions 20 22 35,36]. 

From the theoretical point of view, the B — > Vj decay as well as the low-g 2 region 
of B — > L £ + £~ transitions allow for a systematic inclusion of radiative corrections within 



the QCD factorization approach at leading order in the 1/mb expansion 19 37] . In this 
region, the transition FFs (which still determine a major part of non-perturbative input) 
can be obtained from LCSR estimates alone. As it has been discussed, for instance, in [38], 
the high-g 2 region may also be interesting in order to constrain NP contributions (notably 



7 We have checked that using 90% or 50% correlation, instead — of course — globally changes the very 
value of x 2 , but does not influence the optimal parameter values. A similar comment applies to the number 
of individual LCSR points used in the fit. 
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to the short-distance Wilson coefficients Cg and C\o), and therefore our extrapolations of 
LCSR results for the tensor FFs in that region will be particularly relevant for this purpose. 
In the following subsection, we present the results of fitting the specific FFs to both, the SE 
and SSE parameterisations, using LCSR and Lattice data where appropriate as discussed 



in Section 4.1 For the light meson masses, we use ttik = 494 MeV, ttik* = 892 MeV, 



m p = 776 MeV, m = 1.02 GeV. 

B — >- K form factors: In Figs. [l}|5j we show the fit for the various B — )■ K FFs, which 
enter, for instance, the radiative B — > K£ + £~ and B — > Kvv decays. We compare the 
result of the SE and SSE parameterisations using LCSR data, and investigate the changes 
when the Lattice data is included. The numerical results for the best-fit parameters of the 
SE and SSE fit are found in corresponding Tables |5}|6} The covariance matrices for these 
fits can also be found in Appendix [E| 

Generally, both parameterisations are seen to fit the data well, and importantly, we find 
agreement with the Lattice predictions for Av,o an d Av,t, even when they are not included 
in the fit. We therefore consider our extrapolation of LCSR data for the tensor FF At,o to 
the high-g 2 region, where Lattice data does not exist, as sufficiently reliable. The quality 
of the fits is astonishingly good, considering the \ 2 values for only two free parameters in 
the expansion. The differences between the SE and SSE are only marginal, which can be 



traced back to the usage of the optimised value for the auxiliary parameter t in (16). The 
dispersive bounds turn out to be far from being separated, and therefore they have only 
little impact on the FF fit. This observation is in line with other studies of heavy-to-light 



FFs in the literature, see e.g. 14,18,39 



In order to address the question of potential contributions from higher-order terms in 
the SE, we consider the LCSR prediction for the vector and tensor form factors .4.^0, ^t,o 
and fit to a SE with K = 3, where the coefficient of the z 2 term in the expansion 
has been fixed to values between [—0.9, +0.9], representing almost the maximal range 
allowed by the unitarity constraints. The results of these fits are shown in Figs. 2j6 



Let us first consider the case, where only LCSR data is used in the fit. We observe 
that — as expected — the constraints from the LCSR points at low values of q 2 are not 
sufficient to determine the behaviour at large q 2 , if higher-order terms in the SE (with no 
further phenomenological constraints) are allowed for. However, the behaviour of the form 
factor corresponding to the extreme values of a.2 does not appear very realistic (even a 
rough numerical estimate of coupling constants for the first low-lying resonances with the 
considered hadronic transition would be sufficient to exclude the curves at the margin). 
Therefore the associated error estimate appears too pessimistic to us. Moreover, typically, 
exclusive semi-leptonic branching fractions are suppressed relative to the inclusive ones by 
a factor of about 20. This suggests that the right-hand side of the dispersive bounds should 
not exceed a value of 5% (instead of 100%). Correspondingly, a realistic range of allowed 
<22-values should rather be taken as [—0.25, +0.25], which is indicated by the thick central 
lines in Fig. [2] The associated uncertainty is only slightly larger than the one obtained 
from the variation of (ao,ai) for the SE fit with K = 2. If, on the other hand, the lattice 
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information (in the case of Av,o) is taken into account, the cases with | ctg j > 0.25 actually 
do not yield a satisfactoring fit anymore, and again the error estimate of the linear fit seems 
to be sufficient to estimate the fit errors. In view of the generic difficulties in estimating 
systematic theoretical uncertainties, we thus consider the SE/SSE parameterisation with 
K = 2 and the associated error estimates as sufficiently reliable for practical applications. 
Another comment applies to the scalar FFs Av,t'- As shown in Table [TJ the combined 



heavy-quark/chiral-symmetry limit considered in 26 predicts a scalar B s resonance below 
.B-ff -production threshold (such a state is also favoured by a Lattice computation in |40]). 
On the other hand, the PDG only finds resonances at masses near/above the production 
threshold. We have therefore chosen to compare two variants of the fit, with/without a 
scalar resonance]^] As can be seen, the fit with a scalar resonance from [26] describes the 
combined Lattice/LCSR data significantly better than the fit without a low- lying resonance 
(where in the latter case again the dispersive bounds constrain the FF to lie systematically 
below the Lattice data). However, within the present uncertainties of Lattice and LCSR 
data, this could only be taken as a very indirect argument in favour of a scalar resonance 
in the anticipated mass region. 



Table 5: B — >■ K: Fit of SE parameterisation to LCSR or LCSR/Lattice results, for Av,o 
(X = 1), A v ,t (X = 3) and A T ,o (X = 1). 



A x 


m R 


a 




Fit to 




i 


Av,o 


5.41 


-2.5 x 10- 


-2 


7.2 x 10- 


-2 


LCSR and Lattice 


0.329 




5.77 x 10- 


-3 


A v ,t 




-6.8 x 10" 


-2 


0.20 




LCSR and Lattice 


0.200 




0.129 




A v ,t 


5.72 


-4.5 x 10~ 


-2 


8.9 x 10- 


-2 


LCSR and Lattice 


0.234 




2.99 x 10- 


-2 


Av,o 


5.41 


-2.4 x 10" 


-2 


6.2 x 10- 


-2 


LCSR 


5.07 x 10" 


-3 


4.43 x 10- 


-3 


A v ,t 




-6.7 x 10- 


-2 


0.18 




LCSR 


1.44 x 10- 


-3 


0.111 




A v ,t 


5.72 


-4.8 x 10- 


-2 


0.11 




LCSR 


1.54 x 10- 


-4 


4.34 x 10- 


-2 


At,o 


5.41 


-2.8 x 10- 


-2 


6.0 x 10" 


-2 


LCSR 


3.94 x 10- 


-3 


4.40 x 10- 


-3 



8 Notice that BZ [8] use an effective resonance mass above production threshold to parameterise the 
scalar FFs. 
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Table 6: B — )■ K: Fit of SSE parameterisation to LCSR or LCSR/Lattice results, for Av$ 
{X = 1), A v ,t (X = 3) and A Tfi (X = 1). 



A x 


m R 


a 


«i 


Fit to 




x Cij&i&j 


Av,Q 


5.41 


0.50 


-1.4 


LCSR and Lattice 


0.940 


6.51 x 10" 3 


A v ,t 




0.54 


-1.7 


LCSR and Lattice 


0.904 


0.142 


A v ,t 


5.72 


0.28 


0.35 


LCSR and Lattice 


0.128 


3.15 x 10~ 2 


Av,0 

A v ,t 
A v ,t 
At,o 


5.41 

5.72 
5.41 


0.48 
0.52 
0.30 
0.48 


-1.0 
-1.4 
0.20 
-1.1 


LCSR 
LCSR 
LCSR 
LCSR 


5.15 x 10~ 3 
2.27 x 10" 3 
7.17 x 10~ 5 
8.15 x 10~ 3 


4.04 x 10~ 3 
9.55 x 10~ 2 
5.32 x 10~ 2 
3.06 x 10~ 3 






Figure 1: B — )■ K: Fit of SE (left) and SSE (right) parameterisation to LCSR (top) and 
to LCSR and Lattice (bottom) for Ay,o- The LCSR and Lattice data are shown by black 
points with error bars in the appropriate q 2 range. 
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5 10 15 20 5 10 15 20 



2 2 

q q 

Figure 2: B — > K: Fit of SE parameterisation to LCSR or LCSR/Lattice for Av$ with 
the parameter varied between [—0.9, +0.9] (thin lines). The thick, dark lines show 
a2 = ±0.25. For comparison, on the right we show again the corresponding fit and error 
estimate for a SE truncated after a\. 
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5 10 15 20 5 10 15 20 




5 10 15 20 5 10 15 20 



2 2 

r r 

Figure 3: B K: Fit of SE (left) and SSE (right) parameterisations to LCSR (top) and 
to LCSR and Lattice (bottom) for Ay,t- The LCSR and Lattice data are shown by black 
points with error bars in the appropriate q 2 range. 
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5 10 15 20 5 10 15 20 



Figure 5: B — > K: Fit of SE (left) and SSE (right) parameterisations to LCSR for At,o- 
The LCSR data is shown by black points with error bars. 
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B — >■ p form factors: Our FF fits for B — >■ p transitions, relevant for the radiative 
B — y p7 and B —y pl + l~ decays, are summarized in Figs. [7-13 and Tables [7] and |8j where 
we again compare the fit to SE and SSE parameterisations. As in the case of B — y K 
FFs, we generally observe similar good results for SE and SSE fits, with the dispersive 
bounds again playing only a minor role in restricting the coefficients of the SE/SSE. The 
covariance matrices for the fits can again be found in Appendix [Ej 

Lattice results are restricted to the (axial-)vector FFs, and we again study how the 
fits change when the Lattice data is included: In case of the FF By,o, the uncertainties on 
the Lattice data are rather large, and the fit is in any case dominated by the LCSR points 
at low values of q 2 . Still, we find that the best-fit curve also well describes the central 
values of the Lattice estimates. The situation is somewhat different for By,i, where the 
central values of the Lattice points do not quite agree with the extrapolation of the LCSR 
prediction. The fit is consistent within Lattice uncertainties, but a rather large value of x 2 , 
dominated by the deviations from the Lattice points, is generated. On the other hand, for 
By,i the Lattice data are competitive with the LCSR input, and we can again observe that 
the extrapolation of the LCSR predictions describes the Lattice points very well, while 
inclusion of the Lattice data in this case leads to a very precise FF description. 

In the remaining cases, we again provide the extrapolations for the pseudoscalar and 
tensor FFs from LCSR input, where Lattice data have not been available. Here, it is to 
be mentioned that the uncertainties for the FF £>t,o are quite large, because we had to 
determine the LCSR input values from the difference of two FFs in (12). Of course, it 
would be desirable to directly calculate the FF Bt,o in the LCSR approach which should 
lead to significantly smaller uncertainties for the input data and the extrapolation to large 
values of q 2 . A similar comment applies to the FF By,o- 
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Table 7: B — >■ p: Fit of SE parameterisation to LCSR or LCSR/Lattice results for Bv,o-2 
(X = 1), Bv tt (X = 3) and £ T , - 2 (X = 1). 



B x 


m R 


Po 


Pi 


Fit to 


X&t 


XT. ft 

i 


Bv,o 
B v ,i 

B V ,2 


5.72 
5.33 
5.72 


-8.0 x 10~ 3 
-3.5 x 10~ 2 
-2.5 x 10- 2 


2.5 x 10~ 2 

0.11 
7.8 x 10~ 2 


LCSR and Lattice 


32.1 


1.98 x 10~ 2 


B v ,o 
B v ,i 

B V ,2 


5.72 
5.33 
5.72 


-7.5 x 10- 3 
-3.7 x KT 2 
-2.3 x 10~ 2 


1.4 x 10~ 2 
8.9 x 10~ 2 
5.2 x 10~ 2 


LCSR 


9.56 x lO" 2 


1.28 x 10~ 2 


B v ,t 


5.28 


-3.2 x KT 2 


8.9 x 10~ 2 


LCSR 


3.81 x KT 3 


2.66 x 10" 2 


Bt,o 
B t ,i 

Bt,2 


5.72 
5.33 
5.72 


-1.4 x 10~ 2 
-1.0 x 10~ 2 
-6.3 x 10~ 3 


-8.3 x 10~ 3 
3.4 x 10~ 2 
1.7 x 10~ 2 


LCSR 


4.18 x KT 2 


1.86 x 10~ 3 



Table 8: B — > p: Fit of SSE parameterisation to LCSR or LCSR/Lattice results for Byfl-2 
(X = 1), Bv tt (X = 3) and B Tfi _ 2 (X = 1). 



B x 


m R 


Po 


Pi 


Fit to 


Xfit 


x CijPiPj 


B v ,o 


5.72 


0.26 


0.14 








B v ,i 


5.33 


0.51 


-1.7 


LCSR and Lattice 


33.0 


1.85 x 10~ 2 


B V ,2 


5.72 


0.40 


-0.15 








Bvfl 


5.72 


0.26 


0.50 








B v ,i 


5.33 


0.54 


-1.4 


LCSR 


4.34 x 10~ 2 


1.10 x 10~ 2 


B V ,2 


5.72 


0.37 


0.24 








B v ,t 


5.28 


0.43 


-1.3 


LCSR 


8.49 x 10~ 3 


2.16 x 10~ 2 


Bt,o 


5.72 


0.35 


0.94 








B t ,i 


5.33 


0.52 


-1.5 


LCSR 


3.57 x 10~ 2 


1.79 x 10~ 3 


Bt,2 


5.72 


0.34 


0.31 
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q 2 I 2 

Figure 7: B — >■ p: Fit of SE (left) and SSE (right) parameterisations to LCSR (top) and 
to LCSR and Lattice (bottom) for By,o- The LCSR and Lattice data are shown by black 
points with error bars in the appropriate q 2 range. 
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q 2 I 2 

Figure 8: B — >■ p: Fit of SE (left) and SSE (right) parameterisations to LCSR (top) and 
to LCSR and Lattice (bottom) for By,i- The LCSR and Lattice data are shown by black 
points with error bars in the appropriate q 2 range. 
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Figure 9: B -> p: Fit of SE (left) and SSE 
to LCSR and Lattice (bottom) for By )2 . Th 
points with error bars in the appropriate q 2 




5 10 15 20 




5 10 15 20 



right) parameterisations to LCSR (top) and 
LCSR and Lattice data are shown by black 
ange. 
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Figure 11: B — > p: Fit of SE (left) and SSE (right) parameterisations to LCSR for Bt,o- 
The LCSR data is shown by black points with error bars. 




f-, 




Figure 12: B — > p: Fit of SE (left) and SSE (right) parameterisations to LCSR for Bt,i- 
The LCSR data is shown by black points with error bars. 




f-, 




Figure 13: B — > p: Fit of SE (left) and SSE (right) parameterisations to LCSR for £>x,2- 
The LCSR data is shown by black points with error bars. 
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B — > K* and B s — >■ <fi form factors: The analysis of B — > K* transitions is more 
difficult on the Lattice, as the if*-meson is unstable. Quenched calculations on the Lattice 
have been attempted for the tensor FFs needed in B — > K*j, at q 2 = 0, but we do not 
include these results in our analysis as the other FFs for this decay have not so far been 
calculated. Therefore, we can only fit to the LCSR data, and our numerical results for the 



best-fit parameters of the SE and SSE fit are found in Tables [9] and [10} The covariance 
matrices for the fits can also be found in Appendix |E} 



Table 9: B — > K*: Fit of SE parameterisation to LCSR results for By,o~2 (X = 1), By,t 
(X = 3) and B T ,o-2 (X = 1). 



B x 


m R 


A) 




Fit to 


Xflt 


^EA 2 

i 




5.83 


-9.4 x 10" 


-3 


1.4 x 10" 2 








B v ,i 


5.41 


-5.0 x 10" 


-2 


0.10 


LCSR 


0.149 


1.86 x 10~ 2 


B V ,2 


5.83 


-3.0 x 10" 


-2 


6.8 x 10~ 2 








B v ,t 


5.37 


-4.4 x 10~ 


-2 


0.11 


LCSR 


1.72 x lO" 3 


4.25 x 10" 2 


Bt,o 


5.83 


-1.9 x 10" 


-2 


-1.9 x lO" 2 








Bt,i 


5.41 


-1.4 x 10- 


-2 


4.1 x 10^ 2 


LCSR 


1.88 x 10" 2 


3.14 x 10~ 3 


Bt,2 


5.83 


-8.0 x 10- 


-3 


2.2 x 10~ 2 









As in the case of B — > K*, Lattice QCD predictions for B s — )■ FFs are lacking, and 
we fit to the LCSR data, only. Our numerical results for the best-fit parameters of the SE 
and SSE parameterisations are found in Tables 11 and 12 The covariance matrices for the 
fits can also be found in Appendix |Ej In all cases, we find a good description of the LCSR 
input at low q 2 , and from the experience in B — > K and B — > p transitions we expect the 
extrapolation to high q 2 to be sufficiently reliable. Still, input from Lattice computations 
- if feasible - for B — > K* and B s — > <p transitions at intermediate values of q 2 would be 
highly welcome. 
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Table 10: B ->■ K*\ Fit of SSE parameterisation to LCSR results for B Vfi -2 (X = l), B v ,t 
(X = 3) and B Tfl _ 2 {X = 1). 



5x 


m R 


h 


h 


Fit to 






B Vfi 


5.83 


0.31 


0.74 








Bv,i 


5.41 


0.62 


-1.4 


LCSR 


5.84 x lO" 2 


1.63 x 10~ 2 


B V ,2 


5.83 


0.45 


0.35 








Bv,t 


5.37 


0.49 


-1.4 


LCSR 


4.63 x 10~ 3 


3.62 x 10~ 2 


Bt,o 


5.83 


0.45 


1.4 








B t ,i 


5.41 


0.60 


-1.5 


LCSR 


9.65 x lO" 3 


2.75 x 10~ 3 


Bt,2 


5.83 


0.42 


0.45 









Table 11: B s — > (f>: Fit of SE parameterisation to LCSR results for B v ,o-2 {X — 1), B v ,t 
(X = 3) and B Tfi -2 (X = 1). 



B x 


m R 


Po 


/3i 


Fit to 


xit 


*EA 2 

i 


Bv,o 


5.83 


-5.8 x 10" 


-3 


3.5 x 10~ 3 








B v ,i 


5.41 


-3.4 x 10" 


-2 


9.6 x 10~ 2 


LCSR 


0.124 


1.29 x 10" 2 


B V ,2 


5.83 


-1.8 x 10" 


-2 


4.7 x HT 2 








B v ,t 


5.37 


-3.4 x 10- 


-2 


9.3 x 10~ 2 


LCSR 


8.93 x 10~ 3 


2.96 x 10~ 2 


Bt,o 


5.83 


-1.1 x 10- 


-2 


-1.5 x 10~ 2 








B t ,i 


5.41 


-9.0 x 10- 


-3 


3.5 x 10~ 2 


LCSR 


4.45 x lO" 2 


1.86 x 10~ 3 


Bt,2 


5.83 


-4.6 x 10- 


-3 


1.4 x 10~ 2 
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Table 12: B s — > (f>: Fit of SSE parameterisation to LCSR results for By,o~2 (X — 1), &v,t 
(X = 3) and £ Ti o_2 (X = 1). 





m R 


A> 


h 


Fit to 


Xlt 


x Yl CijPiPj 




5.83 


0.37 


1.1 










5.41 


0.67 


-2.1 


LCSR 


1.44 x 10~ 2 


1.15 x 10~ 2 


B V ,2 


5.83 


0.50 


0.19 








B v ,t 


5.37 


0.61 


-1.8 


LCSR 


1.57 x lO" 2 


2.58 x 10~ 2 


Bt,o 


5.83 


0.51 


1.7 








Bt,i 


5.41 


0.66 


-2.2 


LCSR 


3.49 x 10- 2 


1.66 x 10~ 3 


Bt,2 


5.83 


0.46 


0.26 









5 Discussion and Conclusions 

We have shown that the form factors (FFs) relevant for radiative and semi-leptonic de- 
cays of B and B s mesons into light pseudoscalar or vector mesons can be conveniently 



parameterised as a series expansion (SE) in the variable z(t) (see the definition in (15)). 
With the current accuracy of theoretical estimates from light-cone sum rules (LCSRs) and 
(where available) Lattice QCD, we found that keeping only two terms in the expansion 
and correctly implementing the analytical behaviour due to below-threshold resonances, 
results in a very good description of the FFs over the whole range of momentum transfer 
in the physical decay region. 

The coefficients of the SE are further constrained by dispersive bounds, exploiting the 
crossing symmetry between the physical 5-meson decay and the pair-production of heavy 
and light mesons by the considered decay current. In order to put the discussion for 
the various FFs on a common footing, we found it convenient to use a FF basis where 
the decay /production currents are projected by transverse, longitudinal and time-like po- 
larisation vectors with respect to momentum transfer t. Considering the corresponding 
projections for the current correlators, the constraints take the simple form as indicated in 
( 46p7 ). We stress that for decays into vector mesons the dispersive bounds constrain the 



sum of (squared) coefficients for the three axial-vector FFs, as well as for the three tensor 
FFs. In a simultaneous fit of all FFs, these constraints are thus stronger than those for 
the individual FFs in that sum. 

In order to determine the correct normalization of the SE, given by the profile functions 
<f>(z(t)), we calculate the current correlators using an OPE, including NLO perturbative 
corrections and the leading non-perturbative contributions from quark, gluon and mixed 
condensates. In particular, we provide the NLO results for the tensor-current correlation 
functions, which are relevant for the FFs appearing in radiative and rare semi-leptonic B 
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decays. 

With these theoretical tools at hand, we have performed numerical fits to LCSR 
(Lattice) predictions at low (medium) momentum transfer for all the FFs appearing in 
B K, p, K* and B s — > <fi transitions. We have also investigated a simplified form of 
the SE, where the profile functions (f)(z(t)) are re-expanded in powers of z(t), while the 
dispersive bounds take a somewhat more complicated form. We find that both the stan- 
dard and the simplified SE give a similarly good description of the FF functions. In those 
cases, where Lattice estimates of the FFs is lacking, the SE is used to extrapolate the 
LCSR predictions to the high-g 2 region. Comparing fits with/without using the available 
Lattice data for B — > K and B — >■ p transitions, we judge these extrapolations to be rather 
reliable. Some of our results could be further improved in the future, by addressing some of 
the following issues: The experimental confirmation of a scalar B s resonance below B-K 
threshold, contributing to the scalar B — >■ K FF. Decreasing the uncertainties in Lattice 
predictions for B — > p axial-vector FFs. The calculation of LCSRs directly in the helicity 
basis. The reliable computation of B — > K* and B s — > <fi FFs on the Lattice. 

In conclusion, we have shown that the parameterisation of heavy-to-light FFs as a 
(truncated) SE in z(t) in combination with theoretical estimates from LCSRs and Lattice 
QCD is very useful, not only for the determination of the CKM element \V u b\ from charged 
semi-leptonic B — >• n or B — >• p decays, but also for the description of FFs for radiative 
and semi-leptonic b — > s and b — » d transitions, which will continue to play a major role 
for the indirect search of new physics effects from rare flavour decays. 
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A Kinematics and Polarization Vectors 

In the following, we consider the rest frame of the decaying 5-meson, with the 3-momentum 
of the final-state meson pointing in the z-direction. The polarisation vectors for a (virtual) 
vector state, with 4-momentum g M = (g°, 0, 0, — \q |), are defined as 

4(<Z) = T^= (0,1,^,0), e^q)^^=(\q\,0,0,-q°), 

e'Aq) = -J=9"- (60) 



For the decay of a B- meson at rest into a light meson with mass m L and momentum k, we 
have in particular 

o P m 2 B -m 2 L + q 2 

g u = m B -E = , q\=k = - , 61 

2m r 2m « 
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with A defined in (|6]). We also define the linear combinations 

ef (?) = = (0, 1, 0, 0) , 3(g) = ^M±^M = (0, 0, i, 0) . (62) 



v/2 



In the same way, the polarisation vectors for an on-shell K* meson with momentum 
(E, 0, 0, \k |) are given as 



4(*0 = T^(0,l,±i,0), ^(fc) = — (|fc|,0,0,£). 



(63) 



B FF Properties 

In the following, we summarize a few useful properties of the helicity-based FFs, following 
from the definitions in (5 p|l0|12 ). From the e.o.m. for vanishing momentum transfer, 
q 2 — > 0, one derives 



and 



^v,o(0)=^v,t(0)=/o(0) = / + (0) 



BvM = BvM = ^o(0) = A 3 (0) , 
Bt,M = Bt i2 (0) = >/2Ti(0) = >/2T 2 (0) 



(64) 



(65) 



while the FFs *4r,(b <£>v,i' ^v,2> an d £>t,o vanish like ^/q 2 . Similarly, at the kinematic 
endpoint q 2 = t_ = (m B — mi) 2 , we obtain the relations 



lim 



£y, 2 (g 2 



lim 



?-+t- B Vfi (q 2 ) 9 2 ^t_ -B r ,o(g 2 



V2, 



(66) 



which has been implemented in our FF parameterisation in Section |4.3| 

In the infinite-mass limit mj — > oo, one can project the heavy quark field in the decay 
current onto its large component in HQET, which results in the well-known HQET spin- 
symmetry relations. With our FF conventions, they read 



[m 



B 



q 2 )A 



[m, 



q 2 )A 



v,t 



(67) 



and 



2m B \[q^B T $ = [m B + q 2 ) B v<c 
2m B \/q 2 B T ,i = (m 2 B + q 2 ) B v ,\ 



(m 



B 



(m 



B 



q 2 ) By, 

q 2 ) B v ,2 



t ■ 



2m B y^B Ti2 = (m| + q 2 ) B V:2 + (m| - q 2 ) B 



v.i 



(68) 



Furthermore, in the limit of large recoil energy to the final state meson, one obtains 



additional relations 41 , 42 that follow from the factorisation of soft and collinear QCD 
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dynamics 43-45 . Excluding radiative corrections for q 2 ,m 2 L <C m 2 B , and using the FFs A 
and £>, these relations read: 



A A m B , 

A v ,o - A v ,t - —j=^ A T ,Q , 



(69) 



and 



B 



V,0 



^ By* ^ 



me 



B 



T,2 



(70) 



It is also useful to write the differential decay rates for B — > P£ + £~ and B — > V£ + £~ 
in the naive-factorization approximation in terms of the new FFs. For the decays into 
pseudoscalar mesons, we obtain the relatively simple expression 

-2 {^rn b Cf) \ 2 i a \ r 2 



dT[5 ->■ P£+r 



dg 2 



2AP 



9 Z 



(^l r> o) 2 + 4 AT 2 ^ ' J A T ,o A v ,o 



+ 2AT 2 (|c^ ff | 2 + |C 10 | 2 ) (A v , f 
where the overall normalization is given by 



(71) 



M=W tb v; sl 



G 2 F a 2 



3 • 2 10 vr 5 m| 



A 3/2 



1/2 



For decays into vector mesons, we decompose the doubly differential decay width as 22,46 
d 2 T[B -> 



dg 2 d cos 9 



[(1 + cos 2 9) H T (q 2 ) + 2 cos 9 H A (q 2 ) + 2 (1 - cos 2 9) H L (q : 



(72) 



Here 9 is the angle between the positively charged lepton and the 3-momentum of the 
B° or I? - -meson, and the functions Ht, Ht, Ht are related to the transverse rate, the 
longitudinal rate and the forward-backward asymmetry, respectively. In terms of the new 
FFs, the naive-factorization approximation for the SM contribution results in 

, 2 



H T (q 2 ) = 2Af : 



(2m b C fy 
q 2 



((£ Til ) 2 + (B T , 2 ) 2 ) 



H L (q 2 



H A (q 2 



r9 Re \2m b Cf Cf*] , 
+ 4Af 2 l - — " JL 9 J (B T ,i Bv tl + B T , 2 B Vt2 ) 



+ 2Af 2 {\Cf\ 2 + \C w \ 2 ) {(B VA ) 2 + {B v , 2 f) , 
2A r2 (2m b Cf ) 2 fa ^ , AKr2 Re[2m b CfCf*} 



{B Tfi f + AN' 



Bt,o By,o 



+ 2Ar 2 (|C 9 eff | 2 + |C 10 | 2 ) (Bvjo) 



8 .V' 2 R ( > { Cf* Cva B Vil B v , 2 + ^= Cf C 10 {B T) i B v>2 + B T , 2 B v ,i) 



(73) 

(74) 
(75) 
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Finally, we note the simple FF dependence of the differential decay width and the longitu- 
dinal polarization fraction Fi of the vector meson in the decays B — > Vvv. In the absence 



of right-handed currents, we obtain 47 
dT[B -> Vvv] 



dq 2 



and 



12^ 2 Cl{(B Vfi ) 2 + (B VA ) 2 + (B v , 2 ) 2 ) 
(gy.o) 2 

" (B Vfi ) 2 + (B VA ) 2 + (B v , 2 ) 2 ■ 



(76) 



(77) 



At the kinematic endpoint Bv,i(t~) vanishes, and (66) implies F^it ) = 1/3, as required 
on general grounds. 



C Calculation of Wilson Coefficients Xi 
C.l Perturbative Contribution 



X 







Figure 14: One- and two-loop diagrams contributing to the correlation function. The 
crossed circle indicates the insertion of the corresponding scalar, vector or tensor currents. 
The counter-term diagrams related to the fermion self-energies are not shown. 

In this section, we will briefly sketch the evaluation of the one- and two-loop diagrams 



see Fig. 14) contributing to the perturbative part of the correlation functions. We will 



specify the necessary number n of subtractions for the scalar, vector and tensor correlators, 
and determine the corresponding values of xf{ n ) from the Taylor expansion of the Wilson 
coefficients at q 2 = 0. This leads to a major simplification in the calculation, which allows 
to eliminate external momenta in propagator denominators and to use tensor reduction 
and recursion relations to express the two-loop integrals in terms of two fundamental 
master integrals. Furthermore, we will follow the procedure explained in |48] and absorb 
the IR-sensitive contributions to the Feynman integrals (in the limit m — > 0) into the 
corresponding condensate terms, such that our results have a finite limit when m — > 0. 
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We will find in useful to present the result in terms of the dimensionless variable 



M -m 
M + m 



(7* 



where M and m are the masses of the heavy and light quark in the loop. We further define 
the functions 



fi(v) 
fs(v) 



l-v 2 



atanhff] 



In 

Li 2 



1 - v 
1 +v 

Av 



l + v) 



1 - V 

1 

V 



In 
Li 2 



1 + v 



1 + v 



In 



Av 



(i-vy 



l-v 
2 

,2\ 



4(1 + ^ atanh 2 M 



(79) 



which are manifestly symmetric under exchange of light and heavy quarks (v — > —v), and 
take finite values in the limits v — > {—1,0, 1}. 

We will quote our results for scalar, vector and tensor currents. The expressions for 
currents with opposite parity can be simply obtained by changing v — > 1/v. Our expressions 



for scalar and vector currents coincide with 48 ; the results for the tensor currents are new. 



Scalar Correlator: For the correlator of two scalar currents, we obtain 

t <' 2 1 ' '' "" 

LO 



X S {n = 2) 
X S (n = 2) 



(3 + v 2 )(3v 2 -l) ^ 1 



64vr 2 (M + m) 2 v A 
a s Cp 1 



8vr 2 M 2 



10) 



NLO 



An 64tt 2 (M + m) 2 v 4 

6 (3A (1 - v 2 ) 2 + (3 + v 2 )(3v 2 - 1)) ( f 2 (1 - v 2 ) - A In 
-fl (llw 4 - 50w 2 + 23) + ft (47t» 4 - 126v 2 + 103) 

+4 / 3 v 2 (5v 2 - 1) + 2 (29v 4 + 65v 2 - 40) 



m + M 



— y 



8n 2 M 2 An 



| - 24 In 


~ M~ 


2tt 2 




+ 1T + 
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Vector Correlator: For the different projections of the correlator of two vector currents, 
we obtain 



Xl{n = l) 
Xl{n = 1) 



LO 



(3 + v 2 )(3v 2 -l) J_ 
' 2 ^ 8tt 2 ' 



647T 2 f z 



52) 



a,C, 



1 



NLO 



47T 647T 2 V 2 
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fl (25t» 4 + Uv 2 - 23) + 2/i (I9v 4 - 6v 2 + 23) 
+4/3 (5v 2 - 1) - 23 + Uv 2 + I3v 4 



v-^l 0> s Cf 1 



1 27T 2 

4tt 8tt 2 V 2 + 3 



(83) 



and 

X V T {n = 2) 



-2lv b + 53v 4 + 13V 



3 v-*i 
— ->■ 



LO 



NLO 



512tt 2 (M + m) 2 V 

q; s Cf 1 
4tt 1536tt 2 (M + m) 2 v 4 
.4 



32tt 2 M 2 ' 



(84) 



-fl (803v 6 - 863v 4 - Ihhv 2 - 73) - 2/i (677v 6 - 741u 4 + 279v 2 + 73) 
-4/3 u 2 (I9w 4 - 86w 2 - 5) + 73 + 323w 2 + 755w 4 - 55 lv 
a s CF 3 



'25 2tt^ 
4tt 32vr 2 M 2 V ~6~ + ~3~ 



(85) 



Tensor Correlator: The relevant projection of the tensor current gives rise to 



Xt(" = 3)| 



LO 



-9/i (v 2 - I) 2 (3v 2 + 1) + 4 (-9v 6 + 21v 4 + v 2 + 3) 
256tt 2 (m + M)V 



Xr(n = 3)| 



4tt 2 M 2 



(86) 



NLO 



Att 384tt 2 (M + m)V 
12 (3(u 2 - l) 2 (3v 2 + 1)/! - 3 - V - 2lV + 9v 6 ) 



x / 2 (l-V)-41n 



m + M 
/i 

-/ 2 (766w 6 - 598w 4 - 142w 2 - 218) 
-/1 (l091u 6 - 1137t; 4 + 297v 2 + 325) 

-8/3 v 2 (7v 4 - 26v 2 - 5) + 107 + 69v 2 + 469v 4 - 325u £ 



10 2tt 2 
4tt 4tt 2 M 2 V 3~ + IT + 11 



s^F 



M 



(87) 



C.2 Condensate Contribution to the Correlation Functions 



In this section we provide the expressions for the contributions of the gluon condensate, the 
quark condensate and the mixed quark-gluon condensate to the various current correlators. 
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The contributions to the coefficient of the scalar and vector correlators to all orders in the 



quark mass and lowest order in the coupling constant can already be found in 48 , and 



we reproduce the results given in that paper. We extend this analysis by determining 
the coefficient functions for the tensor correlators. For the quark and the quark-gluon 



condensate, we employ techniques analogous to that given in 48 and closely follow their 



notation. In case of the gluon condensate, we use the plane-wave technique. 



Quark Condensate and Quark— Gluon Condensate: The starting point for calcu- 
lating the coefficient functions to all orders in the quark masses is a closed expression for 
the non-local quark condensate. The position-space expressions for the projection of the 
non-local quark condensates on the local quark condensate °^ and the local mixed 

quark-gluon condensate (^.g s qaFq\)° read 



:q a (0)qp(x)-)- Fq 



(:^:> (0) r 



4m 



D 



[id + m) 



n=0 



-mV/4) n 



n\T(n + D/2) 



x 



(n — + nm 



n=0 



(-mV/4) n 
\p a n\T(n + D/2) 



Here a and (3 indicate the spinor indices. The corresponding projection of the non-local 
mixed quark-gluon condensate reads 



(:Ma(0)^(0)^(i):) 



1 



X 



^-4(D-l)(D-2)m* vM^:> (0) r(f) 

mV/4) n 



((7/A - 7>A) + mcv) (i0 + m)] V - 



/3« n\T(n + D/2) ' 

n=0 



The relevant diagrams for the contribution of the non-local quark condensate and the non- 



local mixed quark-gluon condensate are given in Figs. [151 and C.2, respectively. The 





Figure 15: Diagrams involving the non-local quark condensate, indicated by the two solid 
dots. The crossed circle symbolises the insertion of the currents. 
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Figure 16: Diagrams involving the non-local mixed quark-gluon condensate, indicated by 
the three solid dots. 

evaluation of the diagrams is simplified by the use of the equations of motion, 



m) (■g s q a (0)F lll/ (0)q^x):) 
(p 2 -m 2 ) (:q(0)q(p):) 



qFq 







:g s qaFq:) 



(0) 



qFq 



2(:gg:> 



(0) 



(:q(0)q(p) ■) 



(i>-m) (:q(0)q(p):) 



(90) 
(91) 
(92) 



Gluon Condensate: For the gluon condensate, it is more convenient to use the so-called 



fixed-point gauge technique, which is described in detail in 49 . In the framework of the 
fixed-point gauge, it is possible to derive an expression for 



{<?k\{j/> + m) + (f + m)a KX ] 



(93) 



which is the basic building block for three lowest-order diagrams shown in Fig. 17 



Figure 17: Diagrams involving the gluon condensate. 
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C.2.1 Results 



Quark Condensate: The quark-condensate contribution to the coefficient for the scalar 
correlation function is given by 

(qq) (v + l) 3 



X b (n = 2) 



qq 



ni (99) 

M 5 



(94) 



8(m + M) 5 w 5 

The same expression (up to an overall normalization factor) is obtained in case of the 
longitudinal projection of the vector correlator, Xl( u = 1)| = (M + m) 2 v 2 x S ( n — 2)| 
The transverse projection of the vector correlator leads to 

__(qq)(v + lf 

qq 



\qq 



Xl{n = 2) 



{7v 2 + l) 



64(m + Mfv 5 
Finally, from the relevant tensor correlator we obtain 

(qq) (v + P 3 



(qq) 



xA n = 3) 



(3^ + 1) 



32(m + M) 5 w 5 



M 5 



(gg) 
M 5 ' 



(95) 



(96) 



Gluon Condensate: The expressions for the gluon-condensate contributions to the var- 
ious xf coefficients read as follows: For the scalar correlator we obtain 



X S \n = 2) 



<^ 2 > 



G 2 



96 (m + M)% 6 



-Qv\l - v 2 ) ( / 2 (1 - ^) - 4 In 



15/i(l -w 2 ) 2 - 15- 4w 2 + 27w 4 
m + M" 



/i 



12M 6 ' (97) 
Again, the same expression is obtained for the longitudinal projection of the vector cor- 
relator, x\( n = 1)|gg = (M + m ) 2 v 2 X S ( n — 2)| GG • For the transverse projection of the 
vector correlator, one has 



Xr(n = 2) 



G 2 



384 (m+ M)V 
45 + 115w 2 + 3w 4 - 195v 6 - 5/i (1 - v 2 ) 2 {2hv 2 + 9) 

~m + M 



+u 2 (l - w 2 ) (35w 2 + 41) ( / 2 (1 - v 2 ) - 4 In 

<fG 2 > 



12M 6 ' 

and for the tensor correlator, we get 

& 2 ) 



(98) 



xA n = 3) 



G 2 



384 (m + M)% 6 
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105 + 91v 2 - 17v 4 - 195v 6 - 5/i (1 - v 2 ) 2 (27v 2 + 17) 

~m + M 



+2(1 - v 2 ) (I5w 4 + 2Av 2 + 5) ( j 2 (1 - - 4 In 



24M 6 

Mixed Condensate: For the mixed-condensate contributions, we finally obtain 

X S in = 2% 

and 

3(qGq) 



(99) 



qGq)(l+v) 3<gGg> 
'2 + f ) — > 



xr(n=i)L 



4(m + M) V> 
<gGg)(l + t;) 



2M 7 ' 



(2 + „) "4 1 



\qG q 4(m + M) 5 v 3 v ~ ' ~ y 2M 5 ' 

where — * s a S a i n proportional to x S ( n — ^)|^ G , as we U as 

X £(n = 2)| 

and 



(100) 



(101) 



gGg)(l+-y) / o 9 \ d^i 3(qGq) 

(35v 3 + 59v 2 + Alv + 9) 4 X ^ X , (102) 



X T > = 3)| 



9Q(m + M) 7 v 
(qGq)(l 



M 7 



15w 4 + 28w 3 + 24v 2 + 12u + 5) 



v -p 7(qGq) 
2M 7 



(103) 



D Decomposition of the tensor-current correlator 

Using the projectors 

q^ p-i)g 2 ' 1 ] 

we decompose the correlator of general tensor currents, 

IW = iJd 4 xe iqx (0\T[q 1 (x)a^q 2 (x)q 2 (0)a aP q 1 (0)}\0) : (105) 

into the two Lorentz-invariant functions H TT and H-lt as follows, 

1 3n TT (g 2 ) 

| g^q v q a + gva^m - 9mW ~ 9^Qa ^ 3iI rT (g 2 ) | n ^ ^2^ 
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where 



P£ a Tlftvap — P£ — ^-Lt{Q 2 ) , 

P£ a P^ H-nvap = ^tt(q 2 ) , 

P^PfU^ = 0. (107) 



In this notation, the correlator of the currents 



3a = q°uaq a q (108) 



leads to q 2 Hz/r(g 2 



E Covariance Matrices 

Here we give the covariance matrices as defined in (58) for the parameters corresponding 
to the best-fit parameters in Tables [5] to 12 

B — >• K form factor fit: 

• The fit of B K FFs to LCSR data alone gives the covariances matrices: 



SE SSE 



A 



no res. 

V.t \ _7 87 v 1 n-4 a no v 1 n-3 I I _j 7t v 1 n-2 



1.56 x 1(T 5 


-1.04 x 10- 


-1.04 x 10~ 4 


9.59 x 10" 4 


1.19 x 10~ 4 


-7.87 x 10- 


-7.87 x 10~ 4 


6.98 x 10~ 3 


6.27 x 10" 6 


-2.72 x 10- 


-2.72 x 10~ 5 


2.19 x 10" 4 


2.1 x 10~ 5 


-6.55 x 10- 


-6.55 x 10~ 5 


5.37 x 10" 4 



4.39 x 10- 3 -2.91 x 10- 2 
-2.91 x 10- 2 0.266 



7.17 x 10- 3 -4.75 x 10- 2 
4.75 x 10- 2 0.423 

, , , , 2.61 x 10" 3 -1.08 x 10" 2 

' * o7o^m-5 omvm-4 ) I -1.08 x 10~ 2 8.86 x 10" 2 

7.63 x 10- 4 6.3 x 10~ 4 

6.3 x 10- 4 8.32 x 10- 3 



For the fit of scalar/vector B — > K FFs to LCSR and Lattice data, we obtain the 
covariance matrices: 



Av,o 



SE SSE 

1.48 x 10- 5 -9.81 x 10- 5 \ / 6.26 x 10~ 3 -4.15 x 10~ 2 
-9.81 x 10- 5 8.76 x 10- 4 j V -4.15 x 10~ 2 0.382 



^no res 



4.82 x 10- 5 -2.03 x 10" 4 \ / 3.08 x lO" 3 -1.39 x 10~ 2 
r ' 1 -2.03 x 10- 4 1.6 xlO- 3 J V -1.39 x 10- 2 0.11 

6.21 x 10- 5 -4.11 x 10~ 4 \ / 3.45 x 10~ 3 -2.37 x 10~ 2 



Ay * \ -4.11 x 10" 4 3.75 x 10- 3 j I -2.37 x 10~ 2 0.261 
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B — >• p form factor fit: 



Fitting to LCSR data alone, the covariance matrices for the B — > p FFs are given 
by: 



SE 



Bv,o 
B v ,i 

B V ,2 

B v ,t 
Bt,o 
B t ,i 

Bt,2 



4.15 x 10~ 6 
-3.2 x 10~ 5 

1.57 x 10~ 5 
-1.28 x 10- 4 

6.45 x 1(T 6 
-5.23 x 10~ 5 

1.19 x 10~ 5 
-9.76 x 10- 5 



-3.2 x 10- 5 
7.93 x 10" 4 

-1.28 x 10~ 4 
1.92 x 1(T 3 

-5.23 x 1(T 5 
7.98 x 10~ 4 



-9.76 x 1(T 5 
1.42 x 10~ 3 

1.01 x lO" 5 1.29 x 10~ 4 
1.29 x 10~ 4 1.72 x 10~ 2 

7.45 x 10~ 7 -4.16 x 10~ 6 
-4.16 x 10~ 6 



2.7 x 10~ 7 
-1.41 x 10~ 6 



5.21 x 10~ 5 

-1.41 x 10~ 6 
1.89 x 10~ 5 



SSE 

5.26 x IO- 3 
-4.33 x 10~ 2 

3.29 x 10~ 3 -2.7 
-2.7 x 10~ 2 



4.33 x 10~ 2 
1.57 



x 10~ 2 
0.396 

1.83 x 10~ 3 -1.42 x 10~ 2 
-1.42 x 10- 2 0.274 

2.19 x 10~ 3 -1.81 x 10" 2 

1 01 . , 1 n-2 n otro 



0.258 

X 10~ 2 
15.3 



1.81 x 10~ 2 

6.49 x 10~ 3 9.63 
9.63 x 10- 2 

1.86 x 10~ 3 -9.82 x 10~ 3 

-9.82 x 10~ 3 0.13 

8.09 x 10- 4 -2.15 x 10~ 3 

-2.15 x lO" 3 6.88 x 10~ 2 



For the fit of vector and axial-vector B — )■ p FFs to LCSR and Lattice data, the 
covariance matrices read: 



B v ,o 
B v ,i 

B V ,2 



SE 



2.62 x 10~ 6 


-1.35 x 10~ 5 


-1.35 x 10- 5 


5.35 x 10- 4 


5.72 x 10~ 6 


-3.08 x 10- 5 


-3.08 x 10- 5 


9.15 x 10~ 4 


1.99 x 10~ 6 


-3.02 x 10~ 6 


-3.02 x 10~ 6 


2.26 x 10~ 4 



SSE 

2.86 x 10~ 3 -5.12 x 10~ 3 
-5.12 x 10~ 3 0.796 

1.24 x 10~ 3 -7.07 x lO" 3 
-7.07 x 10- 3 0.193 

5.21 x 10~ 4 1.52 x 10~ 3 

1.52 x 10~ 3 6.4 x 10~ 2 



B — >• K* form factor fit: The covariance matrices for the B — > K* FFs are given by: 
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SE 



SSE 



Bv,o 
B v ,i 

B V ,2 

B v ,t 
Bt,o 

Bt,1 
Bt,2 





Q lO v 1(1" 


-3.12 x 1(T 5 


1.08 x 10- 


1 QzL v in -5 


— 1 . DO A 1U 


-1.63 x 1(T 4 


3.06 x 10- 


i m y in -5 


— O.OZi A 1U 


-8.52 x 1(T 5 


1.59 x 10- 


1.79 x 10~ 5 


-1.53 x 10" 


-1.53 x 1(T 4 


2.75 x 10- 


/ 1.63 x 1(T 5 


3.6 x 10~ 4 


\ 3.6 x 1(T 4 


2.41 x 10" 2 


1.17 x 10" 6 


-6.23 x 10" 


-6.23 x 1(T 6 


9.82 x 10- 


4.04 x lO" 7 


-2.07 x 10" 


-2.07 x 10~ 6 


3.4 x 10~ 5 



4.85 x 10' 3 
-3.26 x 10~ 2 

3. x 10- 3 
-2.54 x lO" 2 

2.42 x 10- 3 
-1.87 x 10~ 2 

2.24 x HT 3 
-1.93 x lO" 2 



-3.26 x lO" 2 
1.85 

-2.54 x lO" 2 
0.467 

-1.87 x lO" 2 
0.456 

-1.93 x lO" 2 
0.34 

9.38 x 10~ 3 0.246 
17.7 

-1.12 x lO" 2 
0.191 

-2.27 x 10~ 3 
0.11 



0.246 

2.26 x 10- 3 
-1.12 x lO" 2 

1.12 x 10~ 3 
-2.27 x 10~ 3 



B s — y 4> form factor fit: The covariance matrices for the B s — y <j) FFs are given 



B v ,o 
Bv,i 

B V ,2 

B v ,t 
Bt,o 

Bt,1 
Bt,2 



SE 



1.16 x lO" 6 
-8.32 x 10" 6 

8.44 x 10~ 6 
-7.79 x 10~ 5 

3.6 x lO" 6 
-3.26 x 10~ 5 

6.61 x 10" 6 
-6.07 x 10~ 5 

7.03 x lO" 6 
1.75 x 10" 4 1.04 

6.39 x 10~ 7 
-3.91 x 10~ 6 

1.63 x lO" 7 



-8.32 x 10~ 6 
3.47 x 10" 4 

-7.79 x 10~ 5 
1.63 x 10" 3 

-3.26 x 10" 5 
7.08 x 10~ 4 

-6.07 x lO" 5 
1.28 x 10~ 3 

1.75 x 10~ 4 \ 
x lO" 2 J 

-3.91 x 10~ 6 
6.67 x 10~ 5 

-8.97 x lO" 7 



SSE 



4.56 x 10~ 3 
-2.81 x lO" 2 

3.37 x 10~ 3 
-3.15 x lO" 2 

2.91 x lO" 3 
-2.38 x lO" 2 

2.05 x 10- 3 
-1.9 x lO" 2 

1.41 x lO" 2 
0.406 



-2.81 x lO" 2 
1.98 

-3.15 x lO" 2 
0.643 

-2.38 x lO" 2 
0.662 

-1.9 x lO" 2 
0.394 



-8.97 x lO" 7 1.69 x 10" 



0.406 
25.2 

3.37 xlO- 3 -1.93x10-2 
-1.93 x lO- 2 0.35 

1.64 x 10- 3 -3.86 x 10- 3 



-3.86 x 10" 



0.187 



47 



References 



[1] M. Artuso et. ai, B, D and K decays, Eur. Phys. J. C57 (2008) 309-492, 
larXiv: 0801. 18331. 



[2] 
[3] 

[4 

[5] 
[6] 
[7 



[9 
[10 

[11 
[12 

[13 
[14 



M. Antonelli et. ai, Flavor Physics in the Quark Sector, arXiv : 0907 . 5386 



J. A. Bailey et. ai, The B — > niu semileptonic form factor from three-flavor lattice 
QCD: A Model-independent determination of \V ub \, Phys. Rev. D79 (2009) 054507, 
|arXiv:0811.3640| . 

QCDSF Collaboration, A. Al-Haydari et. ai, Semileptonic form factors D — > tt, K 
and II - from a fine lattice. Eur. Phys. J. A 13 (2010) 107-120. 
|arXiv:0903.1664| . 

C. Bernard et. ai, Visualization of semileptonic form factors from lattice QCD, 
Phys. Rev. D80 (2009) 034026, |arXiv: 0906 .2498] . 

E. Dalgic et. ai, B Meson Semileptonic Form Factors from Unquenched Lattice 
QCD, Phys. Rev. D73 (2006) 074502, [hep-lat/060102 1|. 

UKQCD Collaboration, K. C. Bowler, J. F. Gill, C. M. Ma ynard, and J. M. F lynn, 
B plv form factors in lattice QCD, JEEP 05 (2004) 035, |hep-lat/0402023] . 

P. Ball and R. Zwicky, New results on B — > 7T, K,r] decay form factors from 
light-cone sum rules, Phys. Rev. D71 (2005) 014015, hep-ph/0406232] . 

P. Ball and R. Zwicky, Bd, s — » p, u, K* , <fi Decay Form Factors from Light-Cone Sum 
Rules Revisited, Phys. Rev. D71 (2005) 014029, [hep-ph/04 12079] . 

P. Colangelo and A. Khodjamirian, QCD sum rules, a modern perspective, 
|hep-ph/0010175[ 

D. Becirevic and A. B. Kaidalov, Comment on the heavy — > light form factors, Phys. 
Lett. B478 (2000) 417-423, |hep-ph/9904490] . 

J. M. Flynn and J. Nieves, \V u b\ from exclusive semileptonic B — > n decays revisited, 
Phys. Rev. D76 (2007) 031302, |arXiv: 0705. 3553] . 

C. Bourrely, I. Caprini, and L. Lellouch, Model-independent description of B — > -nlv 
decays and a determination of \V u b\, Phys. Rev. D79 (2009) 013008, 
|arXiv:0807.2722| . 

M. C. Arnesen, B. Grinstein, I. Z. Rothstein, and I. W. Stewart, A precision model 
independent determination of \V u b\ from B — > nei>, Phys. Rev. Lett. 95 (2005) 
071802, |hep-ph/0504209|. 



48 



C. G. Boyd, B. Grinstein, and R. F. Lebed, Constraints on form-factors for exclusive 
semileptonic heavy to light meson decays, Phys. Rev. Lett. 74 (1995) 4603-4606, 
|hep-ph/9412324] . 

C. G. Boyd and M. J. Savage, Analyticity, shapes of semileptonic form factors, and 
B tt£u, Phys. Rev. D56 (1997) 303-311, |hep-ph/9702300| . 

I. Caprini, L. Lellouch, and M. Neubert, Dispersive bounds on the shape of 



B ->• D^W form factors, Nucl. Phys. B530 (1998) 153-181, hep-ph/9712417|. 



T. Becher and R. J. Hill, Comment on form factor shape and extraction of \V u b\ from 
B ->■ Phys. Lett. B633 (2006) 61-69, |he p-ph/0509090] . 

M. Beneke, T. Feldmann, and D. Seidel, Exclusive radiative and electroweak b — » d 
and b — >■ s penguin decays at NLO, Eur. Phys. J. C41 (2005) 173-188, 
lhep-ph/0412400] . 

C. Bobeth, G. Hiller, and G. Piranishvili, Angular Distributions of B — >■ K 11 
Decays, JEEP 12 (2007) 040, |arXiv:0709.4174| . 

U. Egede, T. Hurth, J. Matias, M. Ramon, and W. Reece, New observables in the 
decay mode B K*°£ + £~, JEEP 11 (2008) 032, |arXiv: 0807. 2589] . 

W. Altmannshofer et. al, Symmetries and Asymmetries of B — )■ K*fi + fi~ . Decays in 



the Standard Model and Beyond, JEEP 01 (2009) 019, [arXiv : 0811 . 1214 



M. Bartsch, M. Beylich, G. Buchalla, and D. N. Gao, Precision Flavour Physics with 
B -> Kvv and B -> K£ + £~, JEEP 11 (2009) 011, [arXiv: 0909 . 1512| . 

R Ball, \V ub \ from UTangles and B -»■ tc£u, Phys. Lett. B644 (2007) 38-44, 



|hep-ph/0611108 



Particle Data Group, C. Amsler et. al, Review of particle physics, Phys. Lett. 
B667 (2008) 1. 

W. A. Bardeen, E. J. Eichten, and C. T. Hill, Chiral Multiplets of Eeavy- Light 
Mesons, Phys. Rev. D68 (2003) 054024, [hep-ph/0305049] . 

M. A. Shifman, A. I. Vainshtein, and V. I. Zakharov, QCD and Resonance Physics. 
Sum Rules, Nucl. Phys. B147 (1979) 385-447. 

M. A. Shifman, A. I. Vainshtein, and V. I. Zakharov, QCD and Resonance Physics: 
Applications, Nucl. Phys. B147 (1979) 448-518. 

V. A. Novikov, M. A. Shifman, A. I. Vainshtein, and V. I. Zakharov, Operator 
expansion in quantum chromodynamics beyond perturbation theory, Nucl. Phys. 
B174 (1980) 378. 



49 



[30] R. J. Hill, The modern description of semileptonic meson form factors, 
|hep-ph/0606023[ 

[31] C. Bernard et. al, Status of the MILC light pseudoscalar meson project, PoS 
LAT2007 (2007) 090, [arXiv: 0710. 1118] . 

[32] R. Williams, C. S. Fischer, and M. R. Pennington, Extracting the qq condensate for 
light quarks beyond the chiral limit in models of QCD, arXiv: 0704. 2296, 

[33] M. V. Polyakov and C. Weiss, Mixed quark-gluon condensate from instantons, Phys. 
Lett. B387 (1996) 841-847, [hep-ph/9607244| . 

[34] G. D'Agostini, Asymmetric Uncertainties: Sources, Treatment and Potential 
Dangers, |physics"7 0403086, 

[35] C. Bobeth, G. Hiller, and G. Piranishvili, CP Asymmetries in bar B — > K*(— > K7r)££ 
and Untagged B s , B s -»■ <j>(-+ K+R-)££ Decays at NLO, JEEP 07 (2008) 106, 
| arXiv: 0805. 25251. 



[36] A. Bharucha and W. Reece, Constraining new physics with B — > K*fi + fi in the 
early LHC era, |arXiv: 1002 .4310| 

[37] M. Beneke, T. Feldmann, and D. Seidel, Systematic approach to exclusive 
B ->■ V£ + £~,V-f decays, Nucl. Phys. B612 (2001) 25-58, |hep-ph/0 106067] . 

[38] B. Grinstein and D. Pirjol, Precise \V u b\ determination from exclusive B decays: 
Controlling the long-distance effects, Phys. Rev. D70 (2004) 114005, 
|hep-ph/0404250] . 

[39] S. Descotes-Genon and A. Le Yaouanc, Parametrisations of the D — > K£v form 

factor and the determination of g, J. Phys. G35 (2008) 115005, |arXiv: 0804. 0203] . 

[40] UKQCD Collaboration, A. M. Green, J. Koponen, C. McNeile, C. Michael, and 
G. Thompson, Excited B mesons from the lattice, Phys. Rev. D69 (2004) 094505, 
|hep-lat/0312007| . 

[41] J. Charles, A. Le Yaouanc, L. Oliver, O. Pene, and J. C. Raynal, Heavy-to-light form 
factors in the heavy mass to large energy limit of QCD, Phys. Rev. D60 (1999) 
014001, |hep-ph/9812358] . 

[42] M. Beneke and T. Feldmann, Symmetry-breaking corrections to heavy-to-light B 

meson form factors at large recoil, Nucl. Phys. B592 (2001) 3-34, |hep-ph/0008255|. 

[43] M. Beneke and T. Feldmann, Spectator interactions and factorization in B — >• niv 
decay, Eur. Phys. J. C33 (2004) s241-s243, [ hep-ph/0308303| . 

[44] B. O. Lange and M. Neubert, Factorization and the soft overlap contribution to 
heavy- to-light form factors, Nucl. Phys. B690 (2004) 249-278, |hep-ph/03il345 . 



50 



[45] M. Beneke and T. Feldmann, Factorization of heavy-to-light form factors in soft- 
collmear effective theory, Nucl. Phys. B685 (2004) 249-296, |hep-ph/0311335] . 



[46] F. Kriiger and J. Matias, Probing new physics via the transverse amplitudes of 
B° -» AT*°(-» K-T T + )e + £- at large recoil, Phys. Rev. D71 (2005) 094009, 
|hep-ph/0502060] . 

[47] W. Altmannshofer, A. J. Buras, D. M. Straub, and M. Wick, New strategies for New 
Physics search in B — > K*uv, B — > Kvv and B — > X s vv decays, JHEP 04 (2009) 
022, |arXiv: 0902 .0160| . 

[48] M. Jamin and M. Miinz, Current correlators to all orders in the quark masses, Z. 
Phys. C60 (1993) 569-578, hep-ph/920820l] . 

[49] L. J. Reinders, H. Rubinstein, and S. Yazaki, Hadron Properties from QCD Sum 
Rules, Phys. Rept. 127 (1985) 1. 



51 



